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Abstract : 

The  Materials  Research  Laboratory  at  The  Pennsylvania  State  University 
fabricates  diphasic  transducer  composites  made  from  polymers  and  piezoceramic. 
It  is  desired  to  measure  the  loss  factor  of  these  materials  over  the  audible 
frequency  range.  After  consideration  of  a  number  of  measuring  systems,  a 
decision  was  made  to  implement  the  Transfer  Function  Method.  This  method 
consists  of  exciting  a  mass-loaded  rod  into  longitudinal  vibrations.  The 
complex  acceleration  ratio  between  its  ends  is  related  to  the  dynamic  modulus 
by  two,  coupled,  transcendental  equations  derived  by  applying  appropriate 
boundary  conditions  to  the  solution  of  the  longitudinal  wave  equation.  In 
theory,  these  equations  can  be  solved  at  any  frequency  by  an  iterative 
procedure.  In  practice,  convergence  problems  reduce  the  valid  solutions  to 
just  those  obtained  at  the  longitudinal  resonance  frequencies.  Improvements 
made  during  this  study,  enable  the  equations  to  be  solved  at  any  frequency 
except  those  were  the  phase  of  that  ratio  is  zero.  Plots  of  loss  factor  and 
storage  modulus  are  presented  for  a  number  of  composites.  Several  high- 
damping  pol)mieric  materials  were  also  tested.  A  comparison  between  the 
present  results  and  those  quoted  in  other  studies  indicates  that  the  improved 
method  produces  reliable  results  more  quickly  and  with  much  less  effort  than 
previously  required. 
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CHAPTER  1 


INTRODUCTION 


1.1  Problem  Definition 

Recent  Improvements  in  the  development  of  diphasic  polymer 
composites  have  lead  to  the  formulation  of  materials  with  improved 
mechanical  and/or  electrical  properties  as  compared  with  single  phase 
materials.  The  Nippon  Electric  Co..  Ltd  (NEC)  has  succeeded  in  developing  a 
ferrite  compound  which  is  structurally  stiff  (compressive  strength:  150  x 
10®  N/m^),  yet  possesses  desirably  high  damping  (100  times  that  of  iron  or 
aluminum)  ID.  The  material  consists  of  ferrite  particles  embedded  in  a 
polyester  resin.  It  is  postulated  that  the  high  specific  gravity  particles 
resonate  with  the  elastic  resin  causing  high  internal  damping. 

In  a  similar  vein,  the  Materials  Research  Laboratory  (MRL)  at  The 
Pennsylvania  State  University,  has  been  fabricating  diphasic  transducer 
materials  made  from  polymers  and  piezoelectric  ceramics.  When  composite 
transducers  are  driven  at  high  frequencies  they  sometimes  exhibit  a  mode  of 
vibration  in  which  the  component  materials  vibrate  out  of  phase.  Internal 
friction  between  the  filler  and  matrix  materials  causes  the  damping. 
Additional  damping,  due  to  electromechanical  energy  conversion  and  ohmic 
dissipation  may  also  be  realized  if  the  piezoceramic  particles  are  polarized. 

A  method  of  measuring  damping  properties  In  the  frequency  range  0  to 
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20  kHz,  for  a  composite  of  limited  size,  was  necessary  for  acoustic  loss 
comparison.  Since  damping  properties  are  a  strong  function  of  frequency  for 
viscoelastic  materials,  the  method  also  had  to  be  capable  of  computing  the 
damping  values  as  a  function  of  frequency. 

After  consideration  of  various  procedures  to  measure  the  dynamic 
modulus  of  elasticity,  a  technique  using  the  acceleration  transfer  ratio 
between  the  two  ends  of  a  mass  loaded  rod  was  selected  [2-5].  A  long  thin 
sample  is  excited  into  longitudinal  vibration  and  the  relative  amplitude  and 
phase  information  recorded  at  frequencies  in  the  audible  range.  The 
amplitude  ratio  and  phase  at  the  longitudinal  resonances  (phase=  ±90*)  is 
related  to  the  elastic  modulus  and  loss  factor  (t\)  by  a  pair  of  coupled, 
nonlinear,  transcendental  equations  derived  from  solving  the  equation  of 
motion  of  a  mass-loaded  rod  excited  into  longitudinal  vibrations.  In  theory, 
the  equations  can  be  solved  at  any  frequency,  providing  amplitude  and  phase 
information  is  known.  In  practice,  convergence  problems  limit  usable  data 
to  only  the  longitudinal  resonances. 

An  improved  computational  algorithm  has  been  developed  which  can 
yield  loss  and  elastic  modulus  values  at  frequencies  other  than  the  ±  90* 
phase  points.  The  method  incorporates  an  improved  "seed"  formulation 
based  on  the  behavior  of  the  accelerometer  transfer  function.  This 
algorithm  facilitates  the  solution  of  the  equations  while  reducing  the 
computation  time. 
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1.2  Discussion  of  Damping 

The  dynamic  characteristics  of  a  structure  can  be  determined  through 
the  characterization  of  the  mass,  stiffness,  and  damping.  The  mass  can  be 
easily  ascertained  using  static  measurements.  The  dynamic  stiffness  is  a 
complex  mechanism  which  almost  always  requires  a  dynamic  test  setup. 
Damping  relates  the  energy  dissipated  to  the  energy  stored  in  the  material. 

Each  dynamic  characteristic  possesses  an  energy  analogy.  Kinetic 
energy  is  associated  with  the  mass,  potential  energy  is  related  to  the 
stiffness,  and  damping  is  a  measure  of  system  energy  losses.  Damping 
energy  is  usually  dissipated  as  heat. 

The  mechanisms  of  damping  are  many.  Interface  friction,  internal 
friction,  fluid  viscosity,  turbulence,  acoustic  radiation,  eddy  currents,  and 
magnetic  hysteresis  are  some  of  the  major  ones  [6].  Internal  and  interface 
friction  constitute  the  most  common  types  of  damping.  Friction 
mechanisms  include:  plastic  slip,  plastic  flow,  dislocation  movements,  and 
inhomogeneous  strain  in  fibrous  materials  [7].  Due  to  the  polymeric  nature 
of  the  materials  tested  in  this  study,  molecular  curling  and  uncurling  of 
polymers  may  also  be  a  major  source  of  damping. 

The  majority  of  the  materials  tested  in  this  study  can  be  classified  as 
viscoelastic.  Viscoelastic  composites  display  Doth  time  and  temperature 
dependent  properties.  All  linear  viscoelastic  materials  possess  a  complex 
modulus  of  elasticity: 

c*=  E’+  iE"  =  E'(l  ♦  in)  =  E'  e’^ 


(1.1) 


where, 
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E*  ’  complex  modulus  of  elasticity 

E‘  =  storage  modulus  of  elasticity 

E"  =  loss  modulus  of  elasticity 
=  loss  factor 

8  =  phase  lag  between  stress  and  strain 

Typical  plots  of  both  the  storage  modulus  and  loss  factor  versus 
frequency  or  temperature  for  a  viscoelastic  material  are  shown  in  Fig.  1.1 
The  loss  factor  (q)  will  be  the  only  measurement  unit  for  damping  used  in 
this  study.  No  standard  nomenclature  for  damping  exists  at  this  time.  This 
is  due  in  part  to  the  great  variety  of  applications  of  damping  for  both 
engineering  purposes  as  well  as  purely  academic  uses.  For  specimens  with 
uniform  lateral  stress  distributions  the  following  damping  relationships 
exist  [7]: 

n=  tan8=  A/TT=  4^/2t7=  2$/(1-$^)'^^=  Af/f„=  \/\=  1/Q  (1.2) 

where, 

A  =  logarithmic  decrement. 

-  specific  damping  capacity. 

5  =  the  damping  ratio 

Af/f^  =  bandwidth  at  the  half  power  points/natural  frequency. 


LOG  SCALE 


INCREASING  FREQUENCY  (CONSTANT  TEMPERATURE) 
DECREASING  TEMPERATURE  (CONSTANT  FREQUENCY) 


Figure  1 . 1  Dependence  of  Modulus  of  Elasticity,  E’,  and  Loss 
Factor,  T\,  on  Temperature  or  Frequency 
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\  =  resonance  amplification  factor  (induced  force/exciting  force). 

Q  =  measure  of  sharpness  of  a  resonance  peak  and  amplification 
produced  by  resonance. 

Along  with  the  problems  of  not  having  a  standard  nomenclature,  there 
are  problems  with  documenting  stress  distributions  and  shape  factors  of 
the  samples  which  are  tested  for  damping.  The  loss  factor  is  not  truly 
independent  of  stress  amplitude.  As  the  stress  amplitude  increases  the 
damping  property  of  the  material  will  also  increase.  If  the  stress  is 
cyclical  the  number  of  fatigue  cycles  will  have  very  little  effect  upon  the 
damping  up  to  a  stress  level  known  as  the  cyclic  stress  sensitivity  limit 
Above  this  value,  stress  at  approximately  80  per  cent  of  the  fatigue  limit, 
the  damping  will  again  increase  [8]. 

Likewise,  the  relative  dimensions  of  a  specimen  can  affect  the  modulus 
of  elasticity,  E^,  and  thus  also  affect  the  loss  factor.  For  rubber-like 
materials  the  resultant  modulus  of  elasticity  is  given  by  [9]; 

1+  BSMe''  (1.3) 

where  E^^  »  resultant  complex  modulus  of  elasticity. 

B  =  numerical  constant  based  on  volume-stress  function 
and  the  stress  distribution.  (B  usually  equals  2.0) 

5  =  shape  factor  (loaded  surface  area/remaining  surface  area) 
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1.3  Complex  Modulus  Measurement  Techniques 

An  abundance  of  techniques  exists  for  measuring  the  damping  properties 
of  materials.  Two  of  the  classic  methods  are  the  Frequency  Response 
Method  and  the  Reverberation  Method  1 10].  The  Frequency  Response  Method 
entails  vibrating  a  cantilevered  specimen  at  its  resonance  frequencies  and 
recording  the  vibration  amplitudes,  preferably  on  a  dB  scale.  By  taking  the 
frequency  bandwidth  between  the  two  -3  dB  points  relative  to  the  peak,  the 
loss  factor  and  storage  modulus  can  be  calculated  from; 

n=Af„/f„  (1.4a) 

E'=  38.24  pL^  f„W  (1.4b) 

where,  p  ■=  sample  density  (kg/m*). 

L  =  unsupported  length  (cm), 
h  =  thickness  in  the  direction  of  vibration  (cm). 

The  Reverberation  Method  uses  a  measure  of  the  vibration  decay  rate  of 
a  specimen.  To  perform  the  test,  the  sample  must  be  excited  into  vibrations 
at  any  resonance  frequency  until  steady  state  conditions  are  achieved.  When 
the  exciting  force  Is  suddenly  removed,  the  resulting  decay  curve  can  be 
used  to  ascertain  the  loss  factor; 

Ti=2.2/T„f„  (1.5) 

where,  Tjq  •  the  time  in  seconds  for  the  vibration  level  to 
decrease  60  dB. 

Nearly  all  methods  for  determining  the  complex  modulus  of  elasticity 
can  be  grouped  into  one  of  the  following  seven  measurement  categories  [1 1], 
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the  shape,  size  and  homogeneity  of  the  sample  will  greatly  determine  which 
method  should  be  used; 

1.  Transverse  Resonance  Mode  Vibration  -  A  small  bar-like  sample  is 

excited  into  transverse  vibration  at  one  end  using  a  non-contact 
transducer  while  a  second  non-contact  transducer  picks  up  the 
displacement  of  the  other  end.  A  plot  of  displacement  versus 
frequency  is  used  to  determine  the  loss  factor  (-3  dB  bandwidth 
method)  and  the  storage  modulus  is  found  using  Eg.  ( 1 .6). 

E‘ =  4.8n2p(L^yhk^^]  (1.6) 

The  coefficient  k^  depends  upon  the  mode  number  and  damping  method. 

2.  Forced  Vibration  -  Used  primarily  for  large,  highly  damped  specimens. 

The  specimen  is  fixed  at  the  two  ends  and  driven  transversely,  at  a 
constant  displacement,  at  its  middle  using  a  sinusoidal  force  F.  The 
equations  for  computing  the  loss  factor  and  storage  modulus  are: 

E'  =  [Fcose/y  ♦  2TT^f\Ap]2LVln^  ( 1 .7) 

Ti  =  ( 1  -  y4T7*f \Ap/F2cos0]  tane  ( 1 .8) 

The  parameter  y  is  the  displacement  of  the  bar's  middle,  I  is  the 
second  moment  of  area  of  the  cross  section  of  area  A,  and  0  is  the 
phase  angle  between  force  and  displacement. 

3.  Flexural  Wave  Method  -  A  probe  microphone  is  used  to  scan  the  length  of 


a  transversely  vibrating  beam,  mounted  on  flexible  supports,  to 
determine  the  resonances  from  the  node  locations.  Using  the  input 
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frequency  and  material  parameters  the  storage  modulus  can  be 
computed  from  Eq.  ( 1 .9); 

E*  =  f„MTT*pALVX  1  (1.9) 

The  constant  Is  dimensionless  and  dependent  upon  the  boundary 
conditions. 

4.  Mechanical  Impedance  Method  -  A  small  specimen  is  driven 

longitudinally  at  a  constant  velocity  using  an  impedance  head.  A  force 
vs.  frequency  plot  is  made  and  the  loss  factor  is  determined  from  the 
-3  dB  bandwidth.  The  storage  modulus  at  the  resonances  (n)  and 
antiresonances  (m)  is  found  from  Eqs.  ( 1 . 1 0)  and  ( 1 . 1 1 ). 

E’ =  pl2L  f/n]VK„  (1.10) 

E*  =  p[2Lf^/m-l/2lVK^  (l.ll) 

The  constants  and  are  correction  factors. 

5.  Transfer  Function  Method  -  A  bar  sample,  with  accelerometers  located  at 

each  end,  is  vibrated  longitudinally.  Using  the  amplitude  and  phase  of 
the  ratio  of  acceleration  at  the  ends  and  the  equations  of  motion 
for  a  mass-loaded  rod  the  damping  and  storage  modulus  values  can  be 
computed,  (also  called  the  Phase  Measurement  Method) 

6.  Thin  Fibers  -  Setup  similar  to  the  Transfer  Function  Method  except  the 

sample  is  loaded  in  tension.  The  loss  factor  is  obtained  from  the  -3 
dB  bandwldths  at  the  resonance  frequencies  and  the  storage  modulus  is 
obtained  using  Eq.  (1.12). 
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E‘ =  4Tl2f^\2p  (,  12) 

7.  Vibration  Decay  Method  -  By  suddenly  removing  the  excitation  force, 
after  steady  state  conditions  are  achieved,  the  loss  factor  can  be 
determined  from  Eq.  ( 1 .5)  and  the  storage  modulus  from  Eq.  ( 1 . 1 2). 

Several  limitations  on  the  size  and  compliance  of  the  samples  available 
for  this  study  reduced  the  method  choice  to  a  couple  of  those  listed  above. 
First,  the  sample  size  had  to  be  kept  to  a  minimum  for  fabrication  purposes. 
The  maximum  convenient  size  for  mixing  and  hot  pressing  was 
approximately  8  cm  by  1  cm^  This  eliminated  the  Forced  Vibration  (2)  and 
Flexural  Wave  methods  (3).  Second,  the  samples  would  include  a  large 
number  of  very  high  damping  materials  so  the  resonance  peaks  would  not 
always  be  apparent.  This  effectively  ruled  out  the  “3  dB  measurement 
schemes,  viz,  the  Transverse  Resonance  Mode  Vibration  ( 1 ),  the  Thin  Fibers 
Method  (6),  and  the  Mechanical  Impedance  Method  (4).  And  last,  the  samples 
were  primarily  too  compliant  for  a  technique  requiring  any  amount  of 
sample  rigidity.  The  Flexural  Wave  Method  (7)  would  not  work  well  with 
extremely  soft  samples. 

The  method  finally  chosen  was  the  Transfer  Function  or  Phase 
Measurement  Method.  The  procedure  consists  of  exciting  a  bar-shaped 
specimen  into  longitudinal  vibration  and  recording  the  complex  ratio  of  the 
acceleration,  or  displacement,  between  the  two  ends.  This  decision  was 
reaffirmed  by  the  findings  of  a  survey  of  current  dynamic  modulus 


n 

techniques  that  was  undertaken  in  ;982  by  Bolt  Beranek  and  Newman  Inc. 
(121. 

That  report  ranks  the  available  systems  for  testing  the  dynamic 
modulus  of  elastomers  according  to  three  cost  ranges.  A  preliminary  list  of 
minimum  requirements  which  all  sytems  had  to  meet  was  described.  Some 
of  the  more  important  of  these  are: 

1.  The  system  measures  at  least  one  frequency  between  200  and  8000  Hz. 

2.  The  accuracy  of  the  system  is  within  ±20%  in  modulus  or  loss  factor. 

3.  The  technique  must  be  capable  of  measuring  two  different  materials. 

4.  The  determination  of  modulus  and  loss  factor  takes  less  than  10  days. 

A  figure  of  merit  was  awarded  to  each  system  based  upon  a  set  of  basic 

attributes.  Typical  attributes  included:  frequency  and  temperature  range, 
accuracy,  time  per  measurement,  availability,  range  of  sample  size,  and 
capability  of  computer  interface.  Table  1.1  assigns  the  attribute  weights 
using  three  schemes.  The  attributes  in  the  “average"  column  are  assigned 
values  according  to  rank  of  importance.  The  most  Important  attribute, 
frequency  and  temperature  range,  received  the  most  points  (16.1  and  15.6) 
while  the  other  attributes  received  decreasing  values  according  to  order  of 
importance.  The  weightings  in  the  column  labeled  "flat"  are  assigned  by 
assuming  that  the  attributes  are  nearly  equal  In  Importance  with  the  more 
important  attributes  still  receiving  a  few  extra  points.  The  "highly  skewed" 
column  gives  a  strong  emphasis  on  frequency  and  temperature  range.  In  all 


12 


Table  1.1  -  Attribute  Weightings  from  the  BBN  study  [12]. 


MrMte.  Average  WeiQhtlna  Envelopes  of  Weightings  Proposed 


Compressive 

Load 

Tensile 

Load 

Compressive 

Load 

Tensile 

Load 

“Flat" 

“Highly 

"Flat" 

“Highly' 

Skewed' 

Skewed' 

Frequency  &  Temperature 

Range 

16.1 

15.6 

10.6 

40.0 

10.6 

38.9 

Static  Modulus  Range 

12.1 

11.7 

10.1 

10.0 

10.1 

9.8 

Accuracy 

11.0 

10.6 

9.6 

10.0 

9.6 

9.8 

Maas.  Tima/Measurament 

10.5 

10.2 

9.0 

10.0 

9.0 

9.8 

Avallablllt//R1sk 

9.0 

8.8 

8.0 

10.0 

7.9 

9.8 

Equipment  Reliability 

8.2 

8n 

8.0 

4.0 

7.9 

3.9 

Total  Meas.  Time/Sample 

7.7 

?.5 

7.45 

4.0 

7.4 

3.9 

Frequency  Sonpling 

6.8 

6.6 

7.45 

4.0 

7.4 

3.9 

Effective  Properties 

5.2 

5.0 

7.45 

4.0 

7.4 

3.9 

Static  Load  &  Pressure 

Range 

5.0 

8.U 

7.45 

2.0 

7.9 

3.9 

Computer  Interface 

4.2 

4.0 

7.45 

1.0 

7.4 

1.2 

Sample  Type 

4.2 

4.0 

7.45 

•  1.0 

7.4 

1.2 
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cases  the  total  point  value  of  all  the  attributes  had  to  add  up  to  100. 

Table  1.2  gives  a  summary  of  the  Figures  of  Merit  for  the  evaluated 
systems.  The  systems  are  grouped  into  three  cost  ranges.  "Type"  refers  to 
the  seven  measurement  categories  enumerated  previously.  Several  of  the 
methods  shown  involve  using  the  Transfer  Function  Method.  Notice  that  the 
top  performer  in  the  least  expensive  group,  the  N5WC  Transfer  Function 
Test,  received  nearly  as  many  points  as  the  top  two  performers  in  the  most 
expensive  group. 

Based  on  this  study  and  the  limitations  previously  mentioned,  a 
decision  was  made  to  Implement  the  Transfer  Function  Method  for  the 
present  study.  The  BBN  study  gives  the  following  information  for  the  NSWC 
Transfer  Method: 

Frequency:  200  to  25,000  Hz. 

Accuracy :  ±5% 

Temperature:  0  to  50’C. 

Static  Modulus:  10®  to  10"  dynes/cm^ 

The  Transfer  Function  Method  looked  particularly  attractive  since  all  of 
the  necessary  electronic  measuring  equipment  already  existed  and  could  be 
linked  to  an  HP  9825B  computer.  The  only  piece  of  equipment  which  needed 
to  be  fabricated  was  a  test  stand  to  hold  the  shaker  and  sample.  Chapter 
Two  details  some  of  the  particulars  of  the  Transfer  Function  Method  and 
also  describes  the  improvements  made  to  the  method  during  this  study. 


Table  1.2  -  Figures  of  Merit  and  Relative  Costs  for  Qualifying 
Systems  from  the  BBN  study  [  1 2]. 
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E.B.  Stringiest  3  70.5  69.4  77.2  72.9  71.7  78.1 

Comparison  Technique  1  69.2  64.4  70.9  69.6  64.5  70.9 

Point  Impedance  4  67.4  64.4  70.9  69.6  64.5  70.9 

Transient  Vibration  Expt.  1  66.1  63.5  59.5  66.6  63.6  62.0 

Tuned  Damper  5  59.3  59.4  68.0  58.1  57.2  67.3 


CHAPTER  2 


THE  TRANSFER  FUNCTION  METHOD 

2.1  Theory  and  Overview 

The  technique  of  obtaining  dynamic  material  parameters  using  the 
acceleration  transfer  ratio  between  the  two  ends  of  a  bar-like  specimen  has 
been  utilized  since  J.  L  Qulmby  first  used  longitudinal  beam  vibrations  to 
measure  viscosity  in  1925 1131  At  present,  the  transfer  ratio  and  material 
parameters  can  be  used  along  with  the  well  known  longitudinal  wave 
propagation  equation  (Eq.  (2.1))  to  determine  E'  and  E"  for  a  bar-like  sample. 

(£•  ♦  lE")(d  Vdx^)=  p(d  Vdt^)  (2. 1 ) 

The  loss  factor, T\,  is  simply  E‘7E’. 

Not  all  of  the  longitudinal  vibration  techniques  require  the 
measurement  of  the  acceleration  transfer  function.  Galkiewicz  and  Karasz 
employed  a  longitudinal  resonance  technique  but  used  the  -3  dB  amplitude 
method  to  obtain  the  logarithmic  decrement  [Ml  Likewise,  Ferguson  uses  a 
similar  method  for  finding  the  loss  using  a  triple-bar  composite  resonance 
technique  [151  D.  I.  6.  Jones  et.  al.  used  a  tuned  damper  technique  to 
simplify  both  the  equations  of  motion  and  the  experimental  procedure  for 
determining  the  dynamic  modulus  of  a  short  rod  [16,171 

The  present  Transfer  Function  Method  evolved  as  the  result  of  two 
theoretical  papers.  Tung-MIng  Lee  combined  the  appropriate  boundary 
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conditions  and  the  longitudinal  wave  propagation  equation  to  model  a 
fixed-free  bar  [18].  Later,  D.  M.  Norris,  Jr.  and  Wun-Chung  Young  extended 
the  solution  to  Include  a  bar  with  an  attached  end  mass  [5].  The  solution  of 
a  bar  with  an  end  mass  is  of  practical  importance  since  it  is  common  to 
have  a  fixture,  or  at  least  an  accelerometer,  attached  to  the  undriven  end  of 
the  bar  or  rod.  Norris  and  Young  used  both  graphical  and  computer  methods 
to  solve  the  resulting  two  equations  representing  the  real  part,  Tp ,  and 
Imaginary  part,  T, ,  of  the  acceleration  ratio,  or  equivalently,  the 

displacement  ratio  of  the  two  ends.  Equations  (2.2)  and  (2.3)  detail  these 
two,  coupled,  transcendental  equations.  The  complete  derivation  can  be 
found  in  Appendix  A. 

Tp  =  coshlCQ]  (cosS  -  RC  sinC)  +  RCQ  cosC  sinhlCO]  (2.2) 

T,  =  sinhlCO]  (sinC  ♦  RC  cosC)  ♦  RCO  sinC  coshlCo]  (2.3) 

where,  R  =  mass  ratio,  M/pAL 

M  =  the  end  mass 
C  -  frequency  parameter,  wL/c 
0  =  tan(6/2) 

c  =  complex  phase  velocity,  (E*/p)'''*sec(8/2) 

This  experimental  procedure  has  been  successfully  utilized  by  many 
people  including  Madigosky  [2,19],  Pritz  [3],  and  Capps  [4].  To  obtain  the 
storage  modulus  and  the  loss  factor  the  solutions  for  C  and  8  should  be 
inserted  into  Eqs.  (2.4)  and  (2.5). 


y\  =  tan(8) 


(2.4) 
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E‘  =  p((jL/C)^  C0S^(8/2)  COS(S)  (2.5) 

The  method  seems  straightforward.  Measure  the  transfer  ratio  between 
the  end  displacements  or  accelerations  and  solve  the  above  equations  using 
a  computer  to  obtain  the  complex  dynamic  modulus.  Unfortunately,  several 
problems  exist  that  limit  both  the  accuracy  and  applicability  of  the 
Transfer  Function  Method. 

2.2  Limitations  of  Present  Transfer  Function  Method 

While  the  Transfer  Function  Method  seems  to  be  gaining  popularity  as  a 
dynamic  modulus  testing  procedure,  it  still  suffers  some  shortcomings  in 
the  form  that  most  practitioners  are  using.  These  problems  can  be  grouped 
into  three  categories:  computer  solution  problems,  data  extension  problems, 
and  lateral  vibration  effects. 

2.2. 1  Computer  Solution  of  the  Coupled  Equations 

The  computer  solution  of  Eqs.  (2.2)  and  (2.3)  can  be  complicated.  The 
most  common  approach  Is  to  use  an  iterative  procedure  such  as  the 
Newton-Raphson  technique  [20]  to  simultaneously  solve  the  two  coupled 
equations.  This  technique  searches  for  the  maxima  or  minima  of  the 
equations  by  locating  “he  zeroes  of  the  derivatives  of  the  equations. 

Although  the  Newton-Raiphson  technique  converges  quickly  to  a  root  of 
the  equations,  there  is  no  guarantee  that  it  will  be  the  correct  root.  It  turns 
out  that  the  number  of  possible  roots  for  a  given  value  of  Tp ,  R,  and  T,  is 
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Infinite;  given  any  value  of  1.  a  value  of  0  can  be  found  which  satisfies  Eqs. 
(2.2)  and  (2.3).  Like  the  solution  to  most  eigenvalue  problems,  the  roots  are 
related  to  the  mode  number. 

An  often  used  approximation  has  been  to  chose  a  seed  based  upon 
neglecting  the  damping  and  solving  the  wave  equation  for  the  frequency 
parameter,  t,  at  the  resonance  frequencies: 

£  =  (2n-l)Tl/2(R+l)  (2.6) 

The  mode  number  is  given  by  “n".  Using  Eq.  (2.6)  to  generate  the  iteration 
starting  point,  or  “seed",  the  solution  of  Eqs.  (2.2)  and  (2.3)  at  the  resonance 
frequencies  can  be  achieved  with  some  degree  of  confidence.  Correct 
solutions,  or  roots,  are  usually  obtained  for  the  first  couple  of  modes.  After 
this,  the  Newton-Raphson  technique  will  often  yield  incorrect  solutions. 
Thus,  the  judgment  of  a  skilled  operator  will  usually  be  necessary  to 
discern  the  good  roots  from  the  bad.  To  further  complicate  matters,  a 
certain  amount  of  symmetry  exists  in  Eqs.  (2.2)  and  (2.3),  that  is,  each 
positive  root  has  a  corresponding  negative  root  which  is  close  in  magnitude. 
This  false  root  can  cause  problems  for  investigators  who  are  not  careful 
about  sign  conventions.  Depending  upon  the  sample  size  and  frequency 
range,  a  typical  test  will  yield  no  more  than  five  or  six  data  points. 

Low  damping  materials  provide  an  additional  problem  for  this 
measuring  system.  Whenever  the  phase  angle  of  the  acceleration  ratio 
approaches  0*,  the  accuracy  of  the  results  becomes  questionable.  The 
solution  of  Eqs.  (2.2)  and  (2.3)  will  result  in  loss  factor  values  very  close  to 
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zero.  The  explanation  for  this  can  be  understood  by  looking  at  Eq.  2.3.  When 
the  phase  angle  equals  zero,  T,  will  also  equal  zero.  The  imaginary  portion 
of  the  acceleration  transfer  ratio  equals  zero  at  any  frequency  if  the 
damping  parameter  is  also  equal  to  zero.  Thus,  Eqs.  (2.2)  and  (2.3)  will  have 
two  solutions,  for  the  same  value  of  i.  tan(S/2),  and  R;  one  with  a  loss 
factor  value  near  zero  and  one  at  the  correct  root.  This  is  the  primary 
reason  why  investigators  have  not  attempted  to  solve  the  equations  at 
nonresonance  frequencies. 

2.2.2  Three  Methods  for  Extending  Results 

There  are  several  ways  to  broaden  the  frequency  range  of  the  results. 
The  simplest  method  Involves  varying  the  weight  of  the  end  mass. 
Increasing  the  end  mass  will  shift  the  resonance  frequencies  to  lower 
values.  Conversely,  a  decrease  in  the  end  mass  will  raise  the  resonance 
frequencies  with  an  absolute  limit  corresponding  to  the  resonance 
frequencies  of  a  free-free  bar.  Using  this  procedure,  the  resonances  can  be 
shifted  to  obtain  results  for  a  wide  range  of  frequencies. 

The  second  method  for  increasing  the  number  of  data  points  is  to  shift 
the  resonances  by  varying  the  length  of  a  sample.  This  will  accomplish  the 
same  effect  as  the  first  method.  A  long  sample  will  naturally  have  a  lower 
primary  resonance  frequency  than  a  short  sample. 

Using  the  present  form  of  the  Transfer  Function  Method,  it  is 
advantageous  to  shift  the  resonances  up  in  frequency.  The  chance  of 
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Obtaining  incorrect  i  and  0  values  increases  greatly  for  higher  modes.  By 
shifting  the  first  few  modes  up  in  frequency,  the  range  of  believable  results 
essentially  shifts  up  in  frequency.  In  addition,  materials  which  are  very 
compliant  will  not  provide  valid  results  at  higher  frequencies.  This  is  due 
to  the  mass-loaded  end  accelerations  being  so  small  that  the  noise  floor  of 
the  test  equipment  is  approached.  Some  materials  reach  this  noise  limit  at 
frequencies  below  4  kHz. 

The  final  means  of  extending  the  frequency  range  of  useful  data  is 
applicable  if  the  material  is  viscoelastic.  A  material  which  is  viscoelastic 
will  exhibit  an  increase  in  strain  as  a  function  of  time  and  temperature 
when  subjected  to  a  constant  stress.  By  recording  damping  and  storage 
modulus  data  at  various  temperatures  and  using  the  time-temperature 
superposition  principle  [21],  the  frequency  range  of  the  data  can  be 
extended.  The  principle,  as  shown  in  Fig.  1.1,  is  based  upon  the  theory  that  a 
decrease  in  temperature  at  a  constant  frequency  causes  the  same  change  in 
dynamic  material  properties  as  an  increase  in  frequency  at  a  constant 
temperature.  The  frequency  shift,  corresponding  to  a  change  from  some 
temperature  T  to  a  reference  temperature  T^^,  is  given  by  the  WLF  [21] 
equation; 

log  a^  -  -c,  (T-T^)  /  (c^  ^  T-T„)  (2.7) 

The  shift  factor  a^  is  equal  to  the  shifted  frequency  divided  by  the 
reference  frequency.  The  constants  c,  and  Cz  are  shift  constants  which  are 
characteristic  of  each  viscoelastic  material  at  a  given  T^^.  If  the  glass 
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transition  temperature  of  the  material  is  used  as  the  reference 
temperature,  T^.  then  c,  and  become  the  universal  constants  (c,  =  1 7.44 
and  C2  =  51.6).  The  WLF  equation  is  primarily  applicable  for  temperatures 
between  and  +  lOO'C.  The  use  of  the  WLF  shift  equation  can 

effectively  extend  the  frequency  range  of  the  data  several  decades  higher 
and  lower  from  those  over  which  valid  data  can  be  obtained  at  a  single 
temperature. 

Most  of  the  Transfer  Function  Method  results  quoted  in  the  literature 
have  been  enhanced  using  the  WLF  shift  equation.  The  procedure  can  be  time 
consuming  and  introduces  additional  errors  due  to  inaccuracies  in  finding 
the  shift  constants  and  applying  the  shift  factor  a^  to  the  experimental 

data.  In  spite  of  this,  the  method  stands  as  a  useful  tool  for  obtaining 
dynamic  material  properties  at  frequencies  outside  the  range  of  test 
limits. 

2.2.3  Lateral  Motion  Effect 

The  wave  equation,  Eq  (2.1 ).  is  based  on  the  assumption  that  wave 
propagation  occurs  only  in  the  axial  directions.  This  approximation  is  only 
valid  if  the  lateral  dimensions  of  the  rod  are  small  compared  with  the 
wavelength.  A  better  equation,  which  takes  into  account  the  effect  of 
lateral  motions,  is  the  Pochhammer-Chree  frequency  equation  for 
longitudinal  waves  in  an  infinitely  long  cylinder  [22].  Using  the  boundary 
condition  that  the  cylindrical  surface  is  traction-free  and  neglecting 


22 

torsional  waves,  the  solution  [23]  to  the  frequency  equation  for  longitudinal 
waves  is: 

(k^-q^)^  [Jo(pa)/J,(pa)]  ♦  4k^pq  [Jo(qa)/J,(qa)]  -  (2p/a)(k^+q^)  =  0  (2.8) 

where  a  =  radius  of  the  cylinder 
p  =  ((u^p/E*)  - 
q  =  [(uVG*)  - 

k  =  complex  wavenumber  (co/c) 

Solving  Eq.  (2.8)  for  its  roots  yields  the  cutoff  frequencies  or 
eigenvalues.  Pritz  [24]  and  Snowdon  [9]  used  a  much  less  complicated 
approximation  to  the  Pochhammer-Chree  equation  developed  by  Love  [25], 

The  latera’  m^iion  is  taken  into  account  to  some  extent  by  the  differential 
equation: 

p  [d^u/dt^  -  (u^Tg)*  d^u/dx^dtn  •  e""  d^u/dx^  (2.9) 

where  u(x)  =  axial  displacement  of  the  bar 

=  complex  Poisson’s  ratio  =  y’(  1  ♦  Iti^) 
r^  -  radius  of  gyration  about  the  x  axis 
For  a  sample  with  a  square  cross-section  of  area  A,  r^*  (A/6)'^^ 

Neglecting  the  lateral  motion  of  the  sample  yields  a  storage  modulus 
and  a  loss  factor  which  deviate  from  the  true  values  of  E’  and  T\  when  the 
lateral  dimensions  become  an  appreciable  portion  of  the  wavelength  at 
higher  frequencies.  For  a  bar  subjected  to  a  time  harmonic  displacement, 

Eq.  (2.9)  becomes; 
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E’(i  ♦  in)  =  £•(!+  \r\)[]  -  (p/E-(i  +  in))w^  y‘^(i  +  inv)^  (2.10) 

The  apparent  storage  modulus  and  apparent  loss  factor  can  be  normalized  by 
the  true  values  of  E‘  and  n: 

E,7E'  "  I  -  (p/E-)  u'  («■  Pj)'  ( 1  -  tiv')  (2.  n ) 

n,/n  *11-  (p/E')  u'  (u'  Pj)'  2(\  /ti)l  /  (  E,'/E')  (2. 1 2) 

Assuming  ny=0,  Eqs.  (2.1 1)  and  (2.12)  can  be  placed  in  a  form  suitable  for 
use  as  a  correction  factor  for  the  apparent  values  (  E^'  and  ng)  obtained  by 
using  the  elementary  theory  (Eq.  2. 1 ). 

E'  »  E;  (1  ♦  4tt^  v‘H  (  1-  nv^))  (2.13) 

n  =  n,  /I  M  4tt2  v'H  r/Xf)^  ( 1  -  nv^  *  2\/n)l  (2. 1 4) 

where,  X^.  =  wavelength  ^  ( 1  /f )(EgVp)’''^ 

The  n  on  the  right  hand  side  of  Eq.  (2.14)  can  oe  approximated  by  n, 
with  little  error.  The  assumption  that  \  is  equal  to  zero  can  be  justified 

by  defining  the  real  and  imaginary  portions  of  the  Poisson’s  ratio  in  terms 
of  the  elastic  and  shear  moduli  and  loss  factors  [9]; 

u*  =  (e’'/26’')  -  1  (2.15) 

This  can  be  separated  into  real  and  imaginary  components: 

v»-u'>iu”  (2.16) 

where 


and 


u*  =  (E72G’)l(1  *y\Q)]-  1 


(2.17) 
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u"  =  (E72G  )[(n  -  \)/i  1  +  V)]  (2. 1 8) 

For  plastic  and  rubberlike  materials  the  elastic  loss  factor,  T\.  is  very  close 
to  the  shear  loss  factor,  hg,  [9]  making  Eq.  (2. 1 8)  approximately  equal  to 

zero.  Since  =  u'Vu',  the  Poisson  loss  factor  will  also  be  close  to  zero. 

Without  the  corrections  in  Eqs.  (2.13)  and  (2.14),  the  elementary  theory 
will  give  storage  modulus  results  at  high  frequencies  which  are  too  low 
while  the  loss  factor  values  will  be  too  high. 

2.3  Improvements  to  the  Transfer  Function  Method 

A  better  understanding  of  the  behavior  of  Eqs.  (2.2)  and  (2.3)  is 
instrumental  in  understanding  why  convergence  problems  occur  when  using 
the  Newton-Raphson  technique  or  any  computer  solution  method.  A 
visualization  of  the  two  equations  as  a  function  of  the  frequency  parameter, 
i,  and  the  damping  parameter,  tan(S/2),  is  obtained  by  making  a  three- 
dimensional  grid  plot.  The  plots  in  Figs.  2.1  and  2.2  give  a  compressed  view 
of  the  real  and  imaginary  components  of  the  acceleration  ratio.  For  ?=  20. 
R=1.0  and  tan(8/2)=  2.0  (corresponding  to  the  highest  points  on  the  curves), 
Tj^  and  T,  are  of  the  order  of  10'^  and  10'®  respectively;  i.e.  very  large 

numbers.  Consequently,  these  transfer  function  amplitudes  are  modified 
using  a  logarithmic  scaling  factor  to  compress  both  the  positive  and 
negative  values  so  that  the  definition  of  the  curve  surfaces  can  readily  be 
discerned.  This  scaling  is  performed  as  follows:  the  logarithm  of  the 
positive  values  greater  than  1  are  plotted,  the  logarithm  of  the  absolute 
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logT^ 


-log  |-T„1 


Figure  2. 1  Three-Dimensional  Grid  Plot  of  the  Real  Transfer 
Function  Amplitude  as  a  Function  of  i  and  tan(S/2) 


for  an  Assumed  Mass  Ratio  R=  I.O. 
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Figure  2.2  Three-Dimensional  Grid  Plot  of  the  Imaginary  Transfer 
Function  Amplitude  as  a  Function  of  i  and  tan(S/2) 
or  an  Assumed  Mass  Ratio  R=  1.0. 
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values  of  the  negative  values  less  than  -1  are  plotted  as  negative  numbers, 
and  the  values  between  1  and  -  I  are  plotted  directly. 

Undulations  in  the  plots  are  analogous  to  the  modes  of  the  mass-loaded 
rod.  Notice  that  the  real  and  imaginary  plots  are  out  of  phase  in  the  sense 
that  a  peak  of  one  curve  nearly  corresponds  to  a  zero  of  the  other  curve.  The 
zeroes  of  the  real  transfer  function  amplitude  indicate  the  locations  of  the 
resonances.  The  zeroes  of  the  imaginary  transfer  function  amplitude 
indicate  the  locations  of  the  antiresonances. 

The  plots  of  the  real  and  imaginary  transfer  function  amplitudes  in  Figs 
2.1  and  2.2  show  only  the  positive  £  and  tan(S/2)  axis  values.  Convergence 
to  negative  values  can  be  achieved,  but  would  result  in  negative  loss  factors 
and  storage  moduli,  neither  of  which  have  any  physical  meaning. 

The  benefit  of  seeing  the  transfer  function  plots  is  that  the  regions  of 
convergence  can  be  accurately  located.  The  trick  is  to  turn  this 
visualization  into  a  solution  algorithm.  Another  variable  complicates  the 
problem.  Equations  (2.2)  and  (2.3)  are  functions  of  three  variables;  £,  0,  and 
R.  The  effect  of  increasing  R  will  be  to  shift  the  undulations  down  the  £ 
axis..  This  corresponds  to  the  shift  in  resonances  when  the  end  mass  is 
changed.  Any  “seed"  generating  algorithm  would  have  to  consider  the 
interrelationship  between  all  three  variables  and  the  mode  number. 

Fortunately,  the  transfer  function  equations  are  very  nearly  Independent 
of  the  damping  parameter,  tan(S/2),  for  a  given  frequency  parameter,  £.  The 
modes  are  almost  uniformly  separated  along  the  C  axis  by  the  value  n.  It  Is 
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only  for  small  values  of  damping  that  the  damping  parameter  becomes  a 
strong  function  of  C. 

Using  the  above  properties,  it  would  not  be  difficult  to  program  a  seed 
algorithm  to  locate  the  approximate  values  of  C  by  assuming  a  value  for 
tan(8/2)  and  solving  for  C  knowing  the  value  of  R.  However,  this  would 
require  too  much  computer  time  and  the  seed  values  would  still  be 
approximate  because  tan(S/2)  is  approximated.  In  addition,  this  same 
iteration  procedure  would  need  to  be  repeated  at  each  frequency. 

A  seed  generation  routine  which  Is  Just  as  accurate  yet  employs  much 
less  computer  time  is  a  two  step  polynomial  curve  fitting  procedure  which 
locates  seed  values  quickly  and  accurately,  usually  within  5%  of  the  true 
solution.  Two  curves  are  necessary.  The  first  polynomial  curve  locates  the 
approximate  ?  value  for  the  first  mode  as  a  function  of  R.  A  plot  of  the 
relationship  between  C  and  R  for  the  first  three  modes  can  be  seen  in  Fig. 
2.3.  The  equation  used  in  this  study  to  approximate  the  curve  for  mode  1  is; 

?,-.002R2-.132R+  1.254  (2.19) 

The  quantity  tan(S/2)  is  assumed  to  be  equal  to  0.1.  Even  if  the  true  value 
of  tan(8/2)  is  0.01,  the  C  values  for  the  first  mode  will  only  differ  by  0.1^ 
from  the  C  values  computed  wit.  t  n(8/2)  equal  to  0.1. 

Once  the  true  value  of  ?,  is  obtained,  it  can  be  used  In  the  second 
polynomial  curve  fitting  procedure  to  approximate  the  higher  mode  seed 
values.  This  polynomial  curve  is  based  upon  the  i  versus  n  (mode  number) 
plot  being  a  straight  line  at  higher  frequencies  with  a  slope  very  near  tt. 
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The  plot  In  Fig.  2.4  shows  that  this  is  indeed  the  case.  The  mass  ratio,  R,  is 
assumed  to  be  equal  to  1.0  for  this  plot.  Different  values  for  R  will  not 
significantly  alter  Fig.  2.4.  Notice  that  the  curve  only  shifts  appreciably  for 
very  high  damping  values,  yet  the  slope  still  remains  the  same  as  that  for 
the  lower  damping  curves.  Thus,  knowing  the  first  mode  value  for  1  and  the 
approximate  slope,  it  should  be  simple  to  find  C  for  the  higher  modes.  The 
simple,  yet  effective,  equation  used  to  approximate  for  higher  modes  is; 

=  (n-l)Ti  +  S,/n  (2.20) 

Again,  the  i  values  are  the  more  important  of  the  two  variables  (1,0)  in 
the  transfer  function  equations  since  the  damping  parameter  can  roughly  be 
considered  independent  of  frequency.  By  applying  a  window  on  the 
convergence  limits  about  this  seed,  the  correct  solution  can  be  assured 
barring  discontinuities  in  the  amplitude  and  phase  data. 

The  key  to  improving  the  solution  method  is  that  once  the  solutions  for 
the  resonance  and  antiresonance  frequencies  are  obtained  accurately,  it  is  a 
simple  matter  to  use  linear  interpolation  between  adjacent  modes  to  find 
solutions  at  any  frequency  providing  the  frequency  lies  between  two  modes 
whose  solutions  have  been  determined.  This  scheme  was  employed  in  the 
current  study  to  obtain  results  across  the  entire  frequency  range  of  interest 
instead  of  only  solving  for  t\  and  E’  at  the  ±90*  transfer  function  phase 
points. 

Another  improvement  made  in  this  study,  which  has  been  overlooked  by 
many  investigators,  is  the  use  of  the  lateral  vibration  correction  as  derived 
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Figure  2.4  Mode  Number  Versus  the  Frequency  Parameter,  I ,  at 


Different  Damping  Levels  for  Mass  Ratio,  R=  1.0. 
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by  Pritz  [24],  Equations  (2. 1 3)  and  (2. 1 4)  were  incorporated  into  the 
solution  algorithm  to  improve  the  accuracy  of  the  results.  The  assumption, 
as  discussed  by  Snowdon  [9],  that  Poisson’s  ratio  can  be  considered  a  real 
value  (T\y=  0)  for  plastics  and  rubbers  was  also  assumed  for  this  study. 


CHAPTER  3 


EXPERIMENTAL  PROCEDURE 


3.1  EQuipment  Setup 

The  Transfer  Function  Method  is  particularly  amenable  to  the  use  of 

% 

accelerometers  for  recording  vibration  levels.  In  most  tests  the 
accelerometer  signals  must  be  calibrated  and  the  vibration  source  must 
have  a  feedback  control  to  monitor  the  vibration  levels.  Neither  is  required 
with  the  Transfer  Function  Method  since  only  the  relative  acceleration 
levels  are  needed.  If  the  two  accelerometers  attached  at  the  ends  of  the 
test  sample  are  well  matched  over  the  frequency  range  of  interest,  their 
output  voltage  ratio  will  be  equal  to  the  acceleration  ratio.  In  addition,  any 
variations  in  the  input  acceleration  level  from  the  shaker  will  not  hamper 
the  test  as  the  response  of  interest  is  the  ratio  between  the  accelerations 
of  the  two  ends  of  the  sample. 

The  equipment  setup  used  In  this  study  is  shown  in  Fig.  3.1.  A  bar-like 
sample  is  harmonically  vibrated  along  its  axis,  at  the  frequency  of  interest, 
by  an  Electrodyne  Shaker.  Two  BBN  501  Accelerometers,  mounted  on  either 
end  of  the  sample,  supply  a  voltage  signal  proportional  to  the  acceleration. 
Each  signal  is  then  amplified  before  feeding  into  a  HP  3570A  Network 
Analyzer.  The  Network  Analyzer  digitizes  and  filters  the  signals  before 
computing  the  relative  amplitude  and  phase  between  the  two  inputs. 
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Figure  3  1  Equipment  Setup  for  Transfer  Function  Method  Tests. 
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The  shaker  is  driven  by  a  HP  3330B  Automatic  Frequency  Synthesizer. 
The  HP  3330B  communicates  to  the  HP  3570A  via  an  interface  bus  (HP  IB). 
This  provides  the  HP  3570A  with  the  input  frequency  on  which  to  perform 
the  bandpass  filtering.  The  filter  bandwidth  used  in  this  study  was  10  Hz. 

The  entire  test  setup  is  controlled  by  a  HP  9825B  computer/ controller 
through  the  same  interface  bus.  The  HP  9825B  also  performs  all  of  the 
numerical  calculations,  stores  all  the  data  on  floppy  disks,  and  controls 
both  a  display  monitor  and  pen  plotter  so  that  plots  of  the  results  can  be 
made.  A  thermal  printer  outputs  the  test  results  in  tabular  form  during  the 
testing.  The  computer,  display,  and  pen  plotter  are  shown  in  Fig.  3.2. 

The  shaker-specimen  fixture,  as  shown  in  Fig.  3.3,  is  a  brass  tripod 
with  neoprene  gaskets  to  reduce  lateral  structure-borne  vibration.  The 
fixture  is  small  enough  to  fit  inside  a  Tenney  Temperature  Chamber  should 
tests  at  various  temperatures  be  desired.  The  sample  is  connected  to  the 
shaker  by  a  small  rectangular  brass  stirrup  which  threads  into  the  shaker. 
One  of  the  accelerometers  is  mounted  inside  this  stirrup  while  the  other 
accelerometer  is  attached  to  a  brass  end  mass.  To  permit  magnetic 
materials  to  be  tested,  such  as  the  NEC  Composite,  brass  was  used  to  make 
many  of  the  test  fixture  components  because  it  is  a  nonmagnetic  metal. 

Two  types  of  epoxy  were  used  in  the  study.  The  accelerometers  were 
each  cemented  to  their  respective  brass  fixtures  with  a  fast-setting 
cyanoacrylate  ester  adhesive  called  Loctite  Super  Bonder  495.  This  type  of 
adhesive  is  similar  to  the  "super  glues"  popular  a  few  years  ago.  The 


Figure  3.2  HP  9825B  Calculator.  HP  1 347A  Display  Monitor, 
and  HP  9872A  Pen  Plotter. 
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sample-Drass  interface  was  sealed  using  a  quick-setting  epoxy  sold  under 
the  trade  name  Epoxi-Patch  by  the  Dexter  Corporation.  The  Epoxi-Patch 
adhesive  hardens  to  an  extremely  brittle  compound  within  5  to  10  minutes. 
These  properties  made  It  Ideal  for  this  study  since  it  did  not  require  a  long 
curing  period  and  the  hardened  epoxy  could  be  chipped  off  the  relatively  soft 
samples  without  tearing  or  otherwise  damaging  them. 

The  Network  Analyzer  and  Frequency  Synthesizer  are  extremely 
accurate.  The  Network  Analyzer  is  accurate  in  amplitude  measurements  to 
0.4  dB  and  accurate  to  0.2*  in  phase  measurements  over  the  frequency  range 
of  the  tests.  Both  of  these  errors  are  far  below  the  error  of  the  rest  of  the 
electronic  equipment  and  the  error  associated  with  poor  mounting  of  either 
the  accelerometers  or  samples. 

The  rest  of  the  equipment  setup  was  tested  to  ensure  that  the 
accelerometer  signals  were  well-matched  and  not  distorted  by  any 
differences  in  the  signal  channels.  First  the  accelerometers  were  mounted 
side-by-side  on  a  fixture  attached  to  the  shaker.  The  response  of  each 
accelerometer  was  recorded  on  a  dB  scale  with  1.0  volt  as  the  reference. 
Figure  3.4a  shows  these  two  signals  on  the  same  plot.  They  match  very  well 
over  the  frequency  range  of  interest  considering  that  they  were  mounted  to 
a  thin  plate  at  slightly  different  locations.  Figure  3.4b  Is  a  comparison  of 
the  same  Input  accelerometer  signal  fed  through  the  two  measurement 
channels.  The  signals  compare  so  well  that  it  is  difficult  to  discern  two 
separate  curves.  This  ensures  that  both  channels  of  the  power  source  and 
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Figure  3.4a  Response  of  Each  Accelerometer  When  Mounted  Side-By-Side 
to  a  Fixture  Attached  to  the  Shaker.  (Reference* 1. 0  volt) 

AMPLITUDE  VS.  FREQUENCY 


Figure  3.4b  Response  of  an  Accelerometer  Signal  Passed  Through 
the  Two  Channels.  (Reference*  1.0  volt) 
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preamplifier  are  well-matched. 

3.2  Materials  Tested 

The  twelve  materials  tested  in  this  study  are  listed  in  Table  3.1  The 
first  seven  consist  of  PZT  (PbZrTiOj)  or  TiOj  powder  embedded  In  a  polymer 
matrix.  These  composites  are  classified  as  0-3  composites.  This  indicates 
that  the  PZT  particles  are  completely  encapsulated  within  the  polymer 
matrix.  The  0  denotes  no  connectivity  in  any  direction  for  the  filler  while 
the  3  Indicates  that  the  polymer  possesses  connectivity  In  all  three 
directions.  The  primary  reason  for  using  PZT  in  the  compounds  instead  of  a 
less  expensive,  dense  powder,  such  as  lead,  is  that  PZT  may  provide 
additional  damping  due  to  electromechanical  losses.  A  comparison  between 
PZT  polymer  composites  which  have  been  polarized  and  the  same  unpolarized 
composites  should  reveal  the  extent  of  the  additional  damping,  if  any, 
caused  by  electromechanical  dissipation. 

One  sample  tested  consisted  of  anatase  (Tiop  in  a  polypropylene 
matrix.  Anatase  has  a  similar  morphology  to  PZT  except  that  It  cannot  be 
polarized.  Loss  factor  differences  between  an  anatase-polymer  composite 
and  a  PZT-polymer  composite  of  the  same  polymer  and  volume  percentage  of 
filler  might  Indicate  that  even  an  unpolarized  piezoelectric  material 
contributes  some  amount  of  electromechanical  damping. 

The  amount  of  PZT  embedded  in  the  polymers  varied  from  0  to  60%  of 
the  volume  fraction.  Some  duplicate  samples,  fabricated  at  different  times, 
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Table  3. 1  -  Materials  Tested  and  their  Manufacturers 


Material  Manufacturer 


PbZrTiOj  (PZT501A) 

TiOj  (Anatase) 
Polyethylene  (Marlex6001) 
Polypropylene 
Poly chi  oroprene 
Spurrs  epoxy 
Eccogel  (Eccogel  Zero) 
Polyurethane  (EN-4) 
Neoprene  (6N  35003) 

NEC  Composite 
SCAB 

Aluminum  (6061) 


Ultrasonic  Powders 
Whittaker  Clark  Daniels 
Phillips  Petroleum 
Phillips  Petroleum 
Polysciences 
Polysciences 
Emerson  &  Gumming 
Conap 

B.  F.  Goodrich 
Nippon  Electric  Company 
B.  F.  Goodrich 
Unknown 
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were  made  to  provide  an  indication  of  the  variability  between  two  batches 
of  the  same  composite. 

The  materials  tested  In  this  study  were  fabricated  at  the  Materials 
Research  Laboratory,  on  the  campus  of  The  Pennsylvania  State  University,  as 
part  of  a  related  research  project  to  characterize  various  properties  of 
ferroelectric/ferromagnetic  polymer  composites.  Two  separate  procedures 
were  used  to  form  the  composites  126]. 

The  polyethylene  and  polypropylene  composites  were  made  using  the 
following  procedure.  First  the  PZT  powder  was  sieved  using  a  No.  325  mesh 
(particle  size  less  than  45  microns).  Next,  the  PZT  and  polymer  were  mixed 
together  using  a  Brabender  mixer  at  40  rpm  for  20  minutes.  After  mixing, 
the  materials  were  hot  pressed  at  pressures  ranging  from  5000  to  10000 
psi.  The  temperature  used  during  this  operation  was  160*C  for  the 
polyethylene  and  195’C  for  the  polypropylene.  The  final  step  was  to  cool 
the  samples  to  room  temperature  while  still  under  pressure. 

A  similar  procedure  was  followed  for  the  Eccogel  and  Spurrs  epoxy 
composites.  The  PZT  powder  was  sieved  using  the  same  No.  325  mesh.  The 
PZT  was  mixed  with  the  polymer  using  either  the  Brabender  or,  in  the  case 
of  the  Eccogel  composites,  by  adding  the  PZT  to  the  liquid  polymer.  The  next 
stage  was  to  pour  the  mixture  into  a  Teflon  mold  and  sonicate  it  to  remove 
air  bubbles.  The  mixture  was  then  cured  at  70’C  for  8  hours. 

Some  of  the  materials  tested  cannot  be  considered  purely  viscoelastic 
since  PZT,  an  elastic  material,  is  used  as  a  filler.  Thus,  the  dynamic 
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modulus  of  the  PZT  composites  should  not  be  as  strong  a  function  of 
frequency  or  temperature  as  that  of  the  pure  polymers. 

The  PZT  was  not  polarized  in  any  of  the  samples  tested.  Samples  will 
be  poled  at  a  later  date  due  to  difficulties  in  polarizing  the  present 
composites.  This  difficulty  was  primarily  due  to  the  low  conductivity  of 
the  polymers.  The  addition  of  carbon  black  to  the  composites  was  being 
considered  at  the  time  of  this  thesis. 

3.3  Dynamic  Modulus  Computer  Program 

The  computer  program  used  to  run  the  test  apparatus  and  solve  Eqs. 

(2.2)  and  (2.3)  for  the  loss  factor  and  storage  modulus  was  written  in  HPL 
(Hewlett  Packard  Language).  HPL  is  a  simple  form  of  BASIC  with  short 
commands  and  single  letter  variables.  It  is  primarily  a  scientific  language 
which  performs  very  well  and  has  an  extensive  error  recovery  system,  but 
does  not  allow  easy  interpretation  due  to  the  lack  of  meaningful  variable 
names  available  in  most  computer  languages.  For  example,  instead  of  using 
"LOSFAC'  to  represent  the  loss  factor,  the  array  "D“  is  used.  For  this  reason 
it  is  felt  that  a  more  instructive  description  of  the  program  can  be  achieved 
by  using  the  flowchart  in  Fig.  3.5.  than  by  including  a  copy  of  the  program 
listing. 

The  entire  program  consists  of  four  sections:  parameter  input,  search 
for  resonances,  solution  of  the  equations,  and  data  storage/retrieval.  The 
first  stage  of  the  program  interactively  loads  the  necessary  material 
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Figure  3  5  Computer  Flowchart  of  Program  Used  to  Compute 
the  Loss  Factor,  T\.  and  Storage  Modulus,  E’ 
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parameters.  After  these  are  loaded,  the  phase  and  amplitude  values, 
measured  at  frequency  increments  of  50  Hz,  are  saved  in  an  array  for  use 
later.  The  primary  resonance  is  also  located  during  this  initial  search  and 
array-loading  operation.  Knowing  the  approximate  frequency  of  the  first 
mode  saves  a  great  deal  of  time  during  the  “search"  stage. 

The  search  for  the  ±  90*  phase  points  is  the  next  step.  The  method  of 
searching  is  iterative.  Phase  values  are  read  until  the  ±  90*  point  is  passed. 
The  frequency  step  is  decreased  and  the  search  begins  at  the  last  frequency 
read  before  passing  the  resonance.  In  this  manner  the  resonance  frequencies 
can  be  ascertained  to  within  0. 1  Hz. 

Once  the  resonance  is  located,  the  phase,  amplitude  and  frequency 
values  are  used  as  input  to  the  solution  algorithm.  Before  beginning  the 
solution  process  the  seed  values  5  and  0  are  determined  using  the  procedure 
discussed  in  Chapter  2.  The  next  step  is  to  use  the  Newton-Raphson  solution 
technique  along  with  Eqs.  (2.2)  and  (2.3)  to  solve  for  the  loss  factor  and 
storage  modulus.  If  this  is  not  the  first  mode  then  the  results  for 
frequencies  between  this  mode  and  the  previous  one  are  determined  by  using 
the  values  for  i  and  0  at  the  two  modes  to  interpolate  the  seed  values  for 
the  frequencies  between  these  two  resonances.  The  data  stored  in  50  Hz 
increments  at  the  very  beginning  are  used  along  with  the  Newton-Raphson 
solution  technique  to  find  t\  and  E’  at  frequencies  between  the  modes. 

This  same  procedure  is  repeated  for  the  remaining  modes  until  the 
frequency  exceeds  20  kHz  or  one  of  the  accelerometer  outputs  reaches  the 


46 


noise  limit  of  the  equipment.  Since  the  driven  end  accelerometer  signal 
does  not  vary  appreciably  with  different  samples,  the  amplitude  ratio 
corresponding  to  the  noise  floor  can  be  estimated.  This  value  is  usually 
near  -70  dB.  The  program  also  contains  an  option  to  extrapolate  seed  values 
past  the  last  mode  and  solve  for  E‘  and  T\  for  frequencies  up  to  the  20  kHz 
limit.  The  upper  frequency  limit  of  20  kHz  was  chosen  based  upon  the 
assumption  that  the  epoxy  bond  could  not  be  considered  truly  rigid  at  higher 
frequencies. 

The  results  are  then  stored  onto  a  floppy  disk  so  that  they  may  later  be 
retrieved,  edited,  or  appended.  Several  different  plotting  routines  are  also 
linked  with  the  program.  The  amplitude  and  phase  data  from  the  test  can  be 
plotted  versus  frequency  to  determine  if  there  are  any  problems  with 
equipment  malfunctions  or  if  the  sample  is  improperly  mounted.  These 
plots  can  also  be  used  to  find  the  frequency  at  which  the  noise  level  of  the 
accelerometer  signals  is  reached.  The  noise  level  is  characterized  by 
random  fluctuations  in  both  the  phase  and  amplitude  plots.  Typical  plots  of 
the  amplitude  and  phase  can  be  seen  In  Figs.  3.6a  and  3.6b.  The  material 
tested  is  polyurethane.  Most  materials  are  characterized  by  only  two  or 
three  modes  in  this  frequency  range  instead  of  the  six  shown  here. 

Several  plotting  routines  exist  for  displaying  the  loss  factor,  the 
storage  modulus,  or  both.  The  particular  plotting  routine  used  contains  an 
option  to  fit  a  least  squares  polynomial  cun/e  through  the  data  if  desired. 

The  dynamic  modulus  testing  procedure  has  been  almost  totally 
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Figure  3.6a  Amplitude  of  Acceleration  Ratio  Between  the  Mass- 

Loaded  End  and  the  Driven  End  of  a  Rod-Like  Specimen. 

phase  vs.  fre:ue.ncy 


Figure  3.6b  Phase  of  Acceleration  Ratio  Between  the  Mass-Loaded 
End  and  the  Driven  End  of  a  Rod-Like  Specimen. 
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automated.  The  input  material  parameters  are  entered  into  the  computer 
and  ten  minutes  later  a  plot  of  the  loss  factor  and  storage  modulus  can  be 
obtained. 


CHAPTER  4 


DISCUSSION  OF  RESULTS 


4.1  Introduction 

The  results  are  presented  in  the  form  of  plots  of  both  the  loss  factor 
and  storage  modulus  versus  frequency.  The  ordinate  axes  are  individually 
scaled  to  make  maximum  use  of  the  plot  area.  The  abscissa  (frequency  axis) 
is  not  scaled  and  thus  some  of  the  plots  only  contain  data  for  the  initial 
portion  of  the  frequency  interval. 

The  samples  were  all  approximately  7  to  8  cm  in  length  and  1  cm^  in 
cross-sectional  area.  The  corrections  for  the  shape  factors,  as  computed 
using  Eq.  ( 1 .3),  were  negligible.  All  tests  were  performed  at  room 
temperature,  24*C.  In  order  to  combine  results  taken  using  different  end 
masses  at  different  times,  the  temperatures  had  to  be  the  same.  The 
ambient  room  temperature  was  chosen  because  it  was  fairly  constant.  In 
addition,  tests  were  attempted  at  temperatures  ranging  from  -20'C  to  45*C 
but  the  available  Tenney  Temperature  Chamber  could  not  provide  a  stable 
temperature  for  the  duration  of  any  test.  The  composites  containing  25%  or 
more  PZT  by  volume  could  not  be  considered  true  viscoelastic  materials  and 
thus  application  of  the  time-temperature  superposition  principle  would 
have  been  inaccurate  for  these  samples.  Note  again  that  in  no  case  were  the 
PZT  particles  polarized,  i.e.  they  function  simply  as  rigid  particles  added  to 
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the  polymer  matrix  material.  The  results  are  grouped  according  to  the 
polymer  matrix  of  the  composites. 

4.2  Eccogel  Composites 

The  composites  composed  of  PZT  and  Eccogel  presented  some 
difficulties  during  the  testing  procedure.  Eccogel  is  extremely  compliant;  it 
has  a  consistency  close  to  gelatin.  It  molds  very  easily,  but  tends  to 
contort  upon  setting.  Results  for  frequencies  higher  than  5  kHz  were 
difficult  to  obtain  for  the  Eccogel  composites  because  the  acceleration 
levels  of  the  mass-loaded  end  approached  the  electrical  noise  limit  of  the 
equipment. 

In  spite  of  these  difficulties,  the  test  results  appear  good.  The  first 
plot,  shown  in  Fig.  4. 1 ,  details  the  loss  factor  and  storage  modulus  results 
of  Eccogel  without  the  addition  of  PZT.  As  might  be  expected,  with  a  very 
compliant  material,  the  damping  is  high.  The  loss  factor  peak, 
corresponding  to  the  region  of  transition  between  the  "rubbery"  stage  and 
the  "glassy"  stage  (see  Fig.  1.1)  of  the  material,  occurs  at  1.8  kHz.  The 
storage  modulus  increases  over  the  entire  frequency  range. 

The  addition  of  PZT  to  the  Eccogel  reduces  the  loss  factor  values  and 
increases  the  storage  modulus.  The  addition  of  a  ceramic,  such  as  PZT,  will 
always  tend  to  increase  the  storage  modulus  of  a  less  rigid  polymer. 

Figures  4.2a  and  4.2b  contain  dynamic  modulus  results  for  identical  Eccogel 
samples  containing  10%  PZT  by  volume.  The  term  "Identical"  is  used  here  to 
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Figure  4.2a  -  Loss  Factor  and  Storage  Modulus  Versus  Frequency 


for  an  Eccogel/ 1 0%  PZT  Sample  (No.  1 ). 
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Figure  4.2b  -  Loss  Factor  and  Storage  Modulus  Versus  Frequency 


for  an  Eccogel/ 10%  PZT  Sample  (No.  2). 
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Indicate  that  the  composites  contain  the  same  proportion  of  polymer  and 
filler.  However,  slight  variations  in  the  fabrication  process  may  cause 
“identical"  samples  to  possess  different  mechanical  properties.  While  the 
general  trends  of  the  loss  factor  curves  are  similar.  Fig.  4.2b  displays 
slightly  higher  damping.  The  loss  factor  peaks  at  approximately  the  same 
frequency  as  the  pure  Eccogel  sample.  The  storage  modulus  curves  for  the 
two  materials  agree  well. 

The  addition  of  more  PZT,  25^  by  volume,  Increases  the  storage  modulus 
by  a  factor  of  5  over  Eccogel  alone.  The  plots  of  the  dyhamic  moduli  of 
Identical  samples  are  shown  In  Figs.  43a  and  43b.  The  loss  factor  curves 
differ  from  those  previously  shown.  Initially  the  loss  factors  decrease  but 
then  increase  without  reaching  a  peak.  This  same  pattern  is  repeated  for 
many  of  the  specimens  containing  other  polymers. 

43  Polyethylene  Composites 

The  composites  made  using  polyethylene  as  the  matrix  material 
generally  have  low  damping  properties  compared  to  the  Eccogel  samples  (an 
order  of  magnitude  lower).  A  plot  of  the  loss  factor  and  storage  modulus  of 
undoped  polyethylene.  Fig.  4.4,  shows  a  gap  where  no  results  are  available. 
These  results  were  discarded  because  they  were  deemed  to  be  unreliable. 
When  the  phase  angle  of  the  transfer  ratio  is  close  to  0*  the  correct 
solution  of  Eqs.  (2.2)  and  (2.3)  becomes  very  difficult  to  obtain.  This  was 
discussed  in  detail  in  section  2.2. 1 .  The  loss  factor  curve  in  Fig.  4.4  most 
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likely  peaks  near  9  kHz  and  then  decreases  In  value  until  it  reaches  the  data 
points  at  12  kHz. 

Figures  4.5a  and  4.5b  contain  the  results  for  identical  specimens  of  259? 
PZT  and  polyethylene.  Again,  the  addition  of  PZT  increases  the  storage 
modulus.  There  is  a  fairly  wide  scatter  in  the  results  for  the  loss  factor  in 
both  plots.  In  the  first  plot  a  peak  appears  to  occur  near  8  kHz. 
Unfortunately,  the  second  plot  contains  no  results  at  this  frequency  to 
either  confirm  or  deny  the  existence  of  a  peak.  The  regions  containing  no 
data  points  signify  frequencies  where  the  phase  of  the  transfer  ratio  is  very 
close  to  0*. 

A  discrepancy  exists  in  the  loss  factor  values  for  frequencies  greater 
than  15  kHz.  In  Fig.  4.5a,  the  loss  factor  increases  for  higher  frequencies 
while  in  Fig.  4.5b  it  remains  stable  over  the  same  region.  This  indicates 
that  either  the  material  or  some  aspect  of  the  test  setup  was  different. 

This  phenomenon  of  increasing  loss  factor  also  occurs  in  the  plots  of  other 
tests  with  different  materials. 

Two  possible  mechanisms  might  account  for  the  rapid  increase  in 
damping  for  the  first  polyethylene  sample  with  25%  PZT  (Fig.  4.5a).  One 
explanation  could  be  that  the  epoxy  holding  the  sample  to  the  test  fixture 
became  ineffective  at  the  higher  frequencies  and  caused  an  increase  in 
friction  at  the  epoxy-specimen  interface.  This  would  explain  why  Fig.  4.5b 
did  not  show  the  same  loss  factor  pattern.  Unfortunately,  a  second  test 
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maj&cY  (Mh) 

Figure  4.5a  -  Loss  Factor  and  Storage  Modulus  Versus  Frequency 


for  a  Polyethylene/25%  PZT  Sample  (No.  1 ). 

LOSS  . 


>  2  4  e  8  la  12  14  16  18  20 


FREQULNCT  (W4i) 

Figure  4.5b  -  Loss  Factor  and  Storage  Modulus  Versus  Frequency 
for  a  Polyethylene/25%  PZT  Sample  (No.  2). 
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using  the  sample  from  Fig.  4.5a  was  run  and  the  test  results  match  the  first 
results  fairly  well;  notice  that  two  separate  storage  modulus  lines  can  be 
seen  on  each  plot.  A  second  possibility  is  that  the  sample  in  Fig.  4.5a  was 
more  porous  than  the  sample,  of  the  same  composition,  shown  in  Fig.  4.5b. 
Variations  in  the  mixing  duration  of  the  polymer  and  filler  could  account  for 
this  difference.  Damping,  in  the  form  of  friction,  is  generated  by  the 
expansion-compression  cycle  of  the  pores.  The  same  effect  occurs  if  the 
material  contains  voids.  This  would  explain  the  increased  damping  at  high 
frequehcies  for  one  sample  as  compared  to  the  relatively  stable  damping  of 
a  second  sample. 

Figure  4.6  contains  the  loss  factor  and  storage  modulus  results  for  a 
60%  PZT  and  polyethylene  composite.  This  was  the  most  dense  sample 
tested  (4.75  g/cm^).  The  storage  modulus  curve  is  somewhat  flat  and 
approximately  twice  the  value  of  the  25%  PZT  composites.  A  sharp  peak  in 
the  loss  factor  occurs  at  approximately  15.6  kHz.  It  is  possible  that 
increasing  the  density  shifts  the  rubbery  frequency  range  to  a  higher  value 
and  causes  a  more  narrow  loss  factor  peak. 

4.4  Polypropylene  Composites 

No  samples  with  PZT  were  available  using  polypropylene  as  the  matrix 
material.  The  two  plots  shown  in  Figs.  4.7a  and  4.7b  contain  the  dynamic 
modulus  results  for  a  pure  polypropylene  specimen  and  a  polypropylene  and 
25%  anatase  specimen.  The  damping  for  the  polypropylene  is  fairly  high  and 
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Figure  4.6  -  Loss  Factor  and  Storage  Modulus  Versus  Frequency 
for  a  Polyethylene/60^  PZT  Sample. 
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Figure  47a  -  Loss  Factor  and  Storage  Modulus  Versus  Frequency 
for  a  Polypropylene  Sample. 
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Figure  4.7b  -  Loss  Factor  and  Storage  Modulus  Versus  Frequency 
for  a  Polypropylene/25%  Ti02  Sample. 
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there  is  a  readily  discernable  loss  factor  peak.  Figure  4.7a  shows  that  the 
loss  factor  and  storage  modulus  for  three  tests  using  three  different  end 
masses  agree  fairly  well  with  each  other. 

A  25%  anatase  and  polypropylene  composite  was  used  for  the  tests 
whose  results  are  shown  in  Fig.  4.7b.  The  anatase  must  also  act  to  stiffen 
the  composite  because  the  loss  factor  decreases  and  the  storage  modulus 
increases  over  the  values  for  polypropylene  alone.  The  loss  factor  initially 
decreases  and  then  shows  a  sharp  increase  at  the  higher  frequencies.  It  is 
possible  that  the  loss  factor  peak  occurs  near  18  or  19  kHz.  Another 
possibility  is  that  the  transition  frequency  (between  the  "rubbery"  and 
■glassy"  regions)  occurs  at  a  much  lower  frequency,  say  around  2  kHz. 
However,  this  would  not  explain  the  sharp  peak  near  18  kHz. 

4.5  Polvchloroorene  Composites 

Two  samples  containing  polychloroprene  were  tested.  One  consisted  of 
only  polychloroprene  and  the  other  contained  25%  PZT  by  volume.  Comparing 
Figs.  4.8a  (no  PZT)  and  4.8b,  it  can  be  seen  that  the  decrease  in  damping  due 
to  the  addition  of  PZT  is  compensated  by  the  interaction  between  the 
polychloroprene  matrix  and  the  PZT  particles.  The  PZT  pa'^ticles  must 
resonate  with  the  matrix  material  causing  additional  friction  losses. 

Notice  that  the  storage  modulus  still  increases  with  the  addition  of  PZT 
particles  by  nearly  a  factor  of  three.  Fig.  4.8a  shows  a  distinct  peak  in  the 
loss  factor  near  6  kHz.  The  loss  factor  curve  in  Fig  4.8b  appears  to  follow 
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the  same  pattern  except  that  It  Is  shifted  higher  In  frequency  by  about  2 
kHz.  This  lends  credence  to  the  idea  that  increasing  the  percentage  of  PZT, 
or  any  dense  filler,  will  shift  the  loss  factor  curve  of  the  polymer 
composite  to  a  higher  frequency. 

4.6  Sourrs  Edoxv  Composites 

Spurrs  epoxy  is  another  low-damping  material.  The  loss  factor,  as 
shown  in  Fig.  4.9,  ranges  between  0.03  and  0.045  for  the  two  separate  tests 
conducted.  The  two  sets  of  test  results  show  a  discrepancy  In  the  loss 
factor  at  frequencies  between  9  and  12  kHz.  This  could  be  a  result  of 
variations  in  the  cure  time  of  the  epoxy  used  to  bond  the  specimen  to  the 
test  fixtures  or  it  might  be  related  to  the  difference  in  the  end  mass  used. 
For  low  damping  materials  the  effect  of  insufficent  cure  time  will  be 
magnified.  The  cure  time  was  kept  fairly  constant  so  this  large  a 
discrepany  Is  probably  a  result  of  the  larger  end  mass  causing  an  increase  in 
the  stress  amplitude  through  the  sample.  This  occurs  with  metals 
subjected  to  cyclical  stress  when  the  fatigue  limit  Is  approached  [7].  It  Is 
not  known  if  the  stress  amplitude  is  great  enough  to  cause  plastic  slip  or 
flow  in  the  polymeric  samples. 

The  same  problem  of  loss  factor  discrepancies  also  occurs  with  the 
results  shown  in  Figs.  4.10a  and  4.10b.  Both  composites  used  for  these 
tests  are  composed  of  25%  PZT  and  Spurrs  epoxy.  For  the  higher 
frequencies,  the  storage  moduli  compare  well,  but  the  loss  factors  do  not. 
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Figure  4.10b  -  Loss  Factor  and  Storage  Modulus  Versus  Frequency 
for  a  Spurrs  Epoxy/25%  PZT  Sample  (No.  2). 
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Again,  the  higher-valued  loss  factor  curve  corresponds  to  the  test  using  the 
larger  end  mass.  The  upper  loss  factor  curve  in  Fig.  4 1  Ob  (two  tests  are 
shown  in  this  figure)  contains  a  sharp  peak  near  12  kHz.  It  looks  like  the 
upper  loss  factor  curve  shown  in  Fig.  4. 10a  might  have  a  peak  near  10  kHz. 
However,  neither  of  the  two  lower  curves  shows  a  peak. 

For  granular  materials  the  loss  factor  increases  markedly  with 
increasing  stress  amplitude,  especially  if  the  acceleration  approaches  the 
acceleration  of  gravity  [6].  Unfortunately,  the  definition  of  "granular"  is 
somewhat  vague.  If  the  PZT -polymer  composites  could  be  considered 
granular,  then  tests  run  with  larger  end  masses  should  show  higher  damping 
levels.  This  indeed  was  the  case  with  the  results  shown  in  Figs.  4. 10a  and 
4.10b. 

4.7  Test  Results  for  Additional  Materials 

The  four  additional  polymeric  materials  tested  provided  good  results. 
The  test  conducted  using  aluminum  yielded  only  a  few  valid  points.  The 
measured  storage  modulus  compares  well  with  values  typically  quoted  for 
the  Young's  Modulus  of  aluminum  (7.0  x  10’“  N/m^).  The  measured  loss  factor 
for  the  first  mode  was  0.002  which  is  much  l^irger  than  the  damping  of  most 
aluminum  alloys  (flgium^  0.0001)  (6].  This  indicates  that  the  inherent 
damping  of  the  test  system  dominates  the  results  for  very  low  values  of 
damping.  The  system  losses  consist  of  epoxy  losses,  brass  fixture  losses, 
and  the  losses  caused  by  the  effect  of  the  accelerometer  cables.  These 
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sources  of  additional  damping  effectively  limit  the  present  test  setup  to 
materials  with  a  loss  factor  no  smaller  than  0.002. 

Figure  4. 1 1  a  Is  a  plot  of  the  storage  modulus  and  loss  factor  for  SCAB. 
50AB  Is  an  acronym  for  a  50und  ABsorbIng  material  marketed  by  B.  F. 
Goodrich  Co..  It  consists  of  aluminum  particles  embedded  in  butyl.  SOAB 
had  the  largest  loss  factor  of  any  material  tested;  it  had  a  peak  loss  factor 
near  1 .3.  The  h  and  E'  curves  for  SOAB  are  reasonably  smooth  and 
continuous. 

A  sample  of  the  Nippon  Electric  Company  (NEC)  composite  discussed  In 
Chapter  1  was  tested.  The  loss  factor  and  storage  modulus  plots  are  shown 
In  Fig.  4.1  lb.  Approximate  values  for  the  loss  factor  and  storage  modulus  of 
this  type  of  composite  [27]  are  h*0.35  and  E-  1.8  x  10’  N/m^  The  present 
test  results  give  slightly  higher  damping  values  and  slightly  lower  storage 
modulus  values. 

The  plot  of  the  dynamic  properties  for  the  neoprene  specimeri  resembles 
the  plots  of  the  Eccogel  composites.  The  loss  factor,  shown  In  Fig.  4.12a, 
peaks  near  3  kHz  with  a  value  of  0.78.  This  result  is  quite  believable.  The 
storage  modulus  Is  low  (1.0  x  10^  N/m^)  like  that  of  the  Eccogel  composites. 

The  final  plot.  Fig.  4.12b,  shows  the  dynamic  properties  of  a 
polyurethane  sample.  The  loss  factor  curve  generally  increases  with 
Increasing  frequency.  The  regions  where  the  loss  factor  tends  to  decrease 
correspond  to  frequencies  near  the  0*  phase  points  of  the  transfer  ratio, 
where  valid  data  could  not  be  obtained.  Both  the  loss  factor  and  storage 
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Figure  4.12a  -  Loss  Factor  and  Storage  Modulus  Versus  Frequenc') 
for  a  Neoprene  Sample. 
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Figure  4.12b  -  Loss  Factor  and  Storage  Modulus  Versus  Frequency 
for  a  Polyurethane  Sample. 
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48  Comparison  of  the  Test  Results 

The  loss  factor  and  storage  modulus  test  results  for  the  various  test 
materials  are  compared  to  each  other  and  to  some  storage  modulus  values 
obtained  using  an  ultrasonic  pulse-echo  technique  at  about  3.0  MHz  [26]  in 
Table  41.  The  materials  are  listed  In  order  of  decreasing  loss  factor.  The 
high-damping  materials  (t\>  0.1 )  are  compared  at  2  kHz  and  the 
lower-damping  materials  are  compared  at  5  kHz.  The  ultrasonic  frequency 
storage  modulus  values  should  be  a  little  larger  than  the  storage  modulus 
values  obtained  in  the  acoustic  range  because  the  storage  modulus  slopes 
slightly  upward  in  the  “glassy"  region.  This  is  indeed  the  case  for  materials 
tested  using  both  methods. 

Additional  comparisons  with  published  results  are  shown  In  Table  42. 

It  is  virtually  Impossible  to  compare  results  of  identical  materials  since 
one  blend  of  a  polymer  will  Inevitably  differ  In  composition  from  another. 
Variations  in  aging,  temperature,  or  frequency  can  also  change  material 
properties  considerably.  The  purpose  of  the  comparisons  is  simply  to 
provide  another  check  against  gross  errors  in  the  results. 
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Table  4.1 


Material 


Eccogel 

SCAB 

Eccogel/ 10*  PZT 
Eccogel/25*  PH 
Neoprene 
NEC  Composite 
Polychloroprene 
Polychloroprene/25*  PZT 
Polyurethane 
Polypropylene 
Polypropylene/25*  TiOj 

Polyethylene 
Spurrs  Epoxy 
Spurrs  Epoxy/25*  PZT 
Polyethy1ene/25*PZT 
Polyethylene/60*PZT 


Comparison  of  Loss  Factor  and  Storage  Modulus 
Results  at  Common  Frequencies. 


Frequency  Loss  Factor  Storage  Modulus 

Acoustic  Ultrasonic*  Pulse/Echo 


kHz 

N/m^ 

N/m2 

2.0 

0.87 

3.78  X  10^ 

2.0 

0.86 

5.89  X  10^ 

2.0 

0.80 

7.27  X  10® 

2.0 

0.67 

3.27  X  10^ 

2.0 

0.51 

1.03  X  10^ 

2.0 

0.25 

9.39  X  10® 

2.0 

0.22 

3.19X  10^ 

2.0 

0.19 

1.22X  10® 

2.0 

0.11 

1.42  X  10® 

5.0 

0.080 

3.96  X  10^ 

5.0 

0.063 

7.80  X  10^ 

5.0 

0.035 

3.07  X  10^ 

5.63  X  10^ 

5.0 

0.029 

3.05  X  10^ 

5.46  X  10^ 

5.0 

0.028 

6.99  X  10^ 

5.0 

0.015 

6.94  X  10^ 

8.25  X  10^ 

5.0 

0.011 

1.75X  I0’° 

2.28 X  lO’O 

-  at  3  MHz 
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Table  42  Comparison  of  Loss  Factor  and  Storage  Modulus  Values 
with  Those  of  Other  Studies. 

Material  Frequency  This  Study  Other  Study  Reference 


(kHz) 

n 

E' 

(N/m^) 

n 

E’ 

(N/m^) 

NEC  Composite 

0.41 

1.0  X  10’ 

0.35 

1.8  X  10’ 

[1] 

Polyurethane 

2.0 

0.11 

1.4x  10® 

0.15 

1.2  X  10® 

141 

Polypropylene 

— 

0.08 

4.0  X  10’ 

0.11 

2.6  X  10’ 

(291 

Aluminum 

__ 

0.002 

7.0  X  10'“ 

0.0001 

7.0  X  10'“ 

161 

CHAPTER  5 


CONCLUSIONS  AND  SUMMARY 


5.1  Conclusions 

The  loss  factor  and  storage  modulus  results  quoted  here  compare 
closely  with  those  of  other  researchers  for  similar  materials.  The  results 
for  high-damping  materials  are  more  characteristic  of  viscoelastic 
materials  than  are  the  results  for  low-damping  materials.  This  is  probably 
a  consequence  of  the  PZT  material  being  an  elastic  material  since  most  of 
the  unusual  results  occur  with  the  addition  of  PZT.  A  good  example  Is  the 
difference  between  Eccogel  with  25%  PZT,  Eccogel  with  10%  PZT,  and  pure 
Eccogel.  While  pure  Eccogel  and  Eccogel  with  10%  PZT  have  a  discernable 
loss  factor  peak,  this  peak  disappears  when  the  percentage  of  PZT  is 
increased  to  25%. 

In  every  case  except  one,  adding  PZT  to  a  polymer  had  the  effect  of 
decreasing  the  loss  factor.  Only  with  polychloroprene  did  the  damping  not 
decrease  appreciably  with  the  addition  of  PZT.  This  probably  is  a  result  of 
the  PZT  mass  particles  resonating  with  the  compliance  of  the  polymer. 
Although  Eccogel  Is  also  a  very  high-damping  polymer.  It  is  not  as  dense  as 
polychloroprene  and  thus  it  may  not  provide  enough  frictional  resistance  to 
the  vibrating  PZT  particles.  This  effect  is  most  likely  responsible  for  the 
high  damping  In  the  NEC  composite  where  ferrite  particles  resonate  with 
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the  polyester  resin  matrix. 

The  test  of  an  aluminum  sample  effectively  set  the  lower  limit  for 
damping  measurements.  Any  material  with  a  loss  factor  less  than  0.002  can 
not  he  tested  accurately  hy  the  measuring  system  used  here  because  the 
inherent  system  losses  will  obscure  the  material  damping.  Since  the 
primary  application  of  the  Transfer  Function  Method  is  to  test  viscoelastic 
or  rubbery  specimens  this  lower  limit  should  not  pose  a  problem.  If  ver7 
low  damping  materials  are  to  be  tested,  either  the  sytem  losses  must  be 
decreased  or  another  dynamic  modulus  test  procedure  must  be  used. 

One  disadvantage  of  the  Transfer  Function  Method  is  that  results  can 
not  be  obtained  for  frequencies  lower  the  first  bar  resonance.  For  some  of 
the  low  damping  composites  this  first  resonance  was  as  high  as  3  kHz.  If 
loss  factor  values  could  be  obtained  for  the  very  low  frequencies  then  many 
of  the  transition  peaks  could  be  discerned.  The  Spurrs  epoxy  samples  appear 
to  have  transition  frequencies  lower  than  the  first  resonance.  The 
supporting  evidence  for  this  is  the  decreasing  loss  factor  values  at  low 
frequencies  shown  in  Figs.  4,9,  4. 1  Oa.  and  4. 1  Ob.  The  loss  factor  peak  for 
these  samples  probably  occurs  before  the  first  resonance. 

Most  of  the  composites  containing  25%  or  more  PZT  are  characterized 
by  relatively  flat  storage  modulus  curves.  Typical  loss  factor  and  storage 
modulus  values  for  this  type  of  PZT  are  ri=  0.002,  and  E'=  8  xlO'°  N/m^  [28]. 
The  PZT  may  act  to  stiffen  the  composite  so  that  the  elastic  modulus 
resembles  that  of  a  ceramic  more  than  that  of  a  viscoelastic  material. 
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Some  of  the  polymers  also  show  very  gradually  Increasing  storage  moduli 
(e.g.,  Fig.  4.9).  The  probable  cause  for  this  behavior  is  that  the  polymer  has 
already  entered  the  “glassy"  region.  The  loss  factor  peak  probably  occurred 
at  a  frequency  below  the  primary  resonance. 

5.2  Summary 

For  applications  at  room  temperature  and  in  the  audible  frequency 
range,  the  Transfer  Function  Method  appears  to  be  a  reasonable  method  for 
measuring  the  storage  modulus  and  loss  factor  of  a  viscoelastic  material. 
The  primary  difficulty  in  using  it  is  that  the  solutions  to  Eqs.  (2.2)  and  (2.3) 
become  suspect  after  the  first  few  modes  if  the  starting  iteration  seed  is 
not  close  to  the  correct  solution.  Equations  (2.2)  and  (2.3)  may  look 
Imposing  but  with  the  aid  of  the  three-dimensional  plots  the  behavior  can  be 
readily  understood.  From  this  point  it  was  a  simple  matter  to  develop  an 
algorithm  guaranteeing  convergence  to  the  correct  solution.  The  algorithm 
developed  during  this  study  to  seed  the  equations  to  the  correct  values  is 
simple,  yet  works  well. 

The  other  improvement  made  to  the  Transfer  Function  Method  was  to 
account  for  the  effect  of  lateral  vibrations  in  the  bar.  This  correction  can 
easily  change  both  the  loss  factor  and  storage  modulus  by  10%  or  more  at 
the  higher  frequencies.  It  should  not  be  neglected  as  commonly  practiced 
when  using  the  Transfer  Function  Method. 

Although  the  loss  factor  and  storage  modulus  plots  may  still  be  subject 
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to  some  experimental  error,  they  are  useful  from  the  standpoint  that  a 
comparison  of  the  damping  between  different  polymer  composites  can  be 
made  with  confidence.  The  test  results  also  seem  to  indicate  that  the 
method  is  capable  of  providing  reproducible  results.  These  two 
requirements  were  considered  the  most  important  to  meet  when  the 
decision  to  develop  a  test  procedure  for  measuring  the  damping  properties 
of  the  PZT-polymer  composites  was  made. 

The  primary  attributes  of  the  Transfer  Function  Method,  with  the 
Improvements  made  during  this  study,  are  that  it  is  very  fast,  easy  to 
perform,  can  be  implemented  at  modest  cost,  is  capable  of  testing 
viscoelastic  materials  in  a  wide  range  of  sizes,  and  requires  almost  no 
operator  Judgment. 

5.3  Recommendations  for  Further  Study 

The  test  setup  in  its  present  form  is  useful  but  with  the  proper 
modifications  it  could  become  much  more  versatile  and  accurate.  The  first 
Improvement  that  could  be  undertaken  would  be  to  decrease  the  inherent 
system  damping.  This  could  be  accomplished  by  using  noncontact 
transducers  for  determining  the  acceleration  ratio  of  the  mass-loaded  and 
driven  ends.  If  the  system  losses  could  be  reduced  to  the  losses  associated 
with  air  friction  then  this  test  arrangement  could  be  used  to  test  very  low 
damping  materials  such  as  metals. 
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Another  Improvement  would  be  to  solve  Eqs.  (2.2)  and  (2.3)  for  the  loss 
factor  and  storage  modulus  at  frequencies  below  the  primary  resonance. 

The  same  extrapolation  procedure  used  to  approximate  seed  values  beyond 
the  last  mode  could  also  be  used  to  seed  the  Newton-Raphson  iteration 
procedure  for  frequencies  smaller  than  the  first  mode  resonance  frequency 
provided  the  amplitude  and  phase  data  were  accurate.  This  would  provide 
valuable  low  frequency  data.  The  behavior  of  the  dynamic  modulus  at  low 
frequencies  seems  to  be  Important  because  it  appears  that  for  many 
composites  the  transition  frequency  Is  below  2  kHz. 

One  final  improvement  to  the  method  would  be  to  extend  the  testing 
procedure  to  include  torsional  vibrations.  The  equations  of  motion  are  the 
same  as  for  the  longitudinal  vibrations  except  that  the  material  and 
displacement  parameters  change.  If  both  the  complex  elastic  and  shear 
moduli  could  be  determined,  the  complex  Poisson’s  ratio  could  also  be  found 
by  using  relationships  among  the  various  elastic  constants.  This  would 
probably  necessitate  that  physical  changes  be  made  to  the  testing 
apparatus.  However,  the  software  changes  would  probably  be  minimal. 
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Derivotion  of  equations  used  to  solve  for  loss  factor 
and  storage  modulus. 


A=  cross-sectional  area  of  the  bar. 
M=  the  end  mass. 

E*=  complex  elastic  modulus. 

I<=  complex  propagation  constant. 

F=  force  on  the  cross-section. 

Assume  a  solution  of  the  one  dimensional,  reduced,  wave 
equation  of  the  form; 

u(x)=  C,  coskx  ♦  C2  sinkx  u(x)=  longitudinal  displacement 

of  the  bar. 

u(0)=  Uq  Uo=  displacement  at  x=0  and  t=0. 

The  axial  force  in  the  bar  at  any  position  x  is  given  by. 

F=  -AE*  au/ax 

The  boundary  conditions  are. 
x=0  u(0,t)=  C,=  Uq 

@  x=L  F=  M  a^u/at^=  -Mq^u(L) 

Therefore: 

-AE*  au/ax|^_^= -flcj^u(L)  (A  4) 

Applying  the  solution  for  u(x)  yields. 

AE*kl-Uo  sinkL  ♦  C2  coskL]=  Uq  coskL  ♦  C2  sinkL)  (A  5) 

or 

C2lriu^  SinkL  -  AE*k  coskLl  =  -(AE*kUo  sinkL  ♦  coskL)  (A  6) 


(A.2) 

(A.3) 


thus. 
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Uol  AE*k  sinkL  ♦  coskL] 
AE*k  coskL  -  SinkL 


Uq!  tankL  ♦  7] 
n  -  TlankLl 


.  -v  M 10^ 
where  7= - 

AE*k 


(A. 7) 

(A.0) 


The  ratio  of  displacement  of  driven  end  to  that  of  the  mass 
loaded  end  is  therefore: 


u(o)  Uo  - ya - - 

=  —7T  =  Uni  tankL  ♦  ^1 

u(L)  u(L)  Uq  coskL  ♦  -  sinkL 

[1  -  TtankLl 


(A.9) 


1-  I'tankL 

[(  coskL  -  VsinkL)  ♦  (  tankL  sinkL  'J'sinkDl 


(A. 10) 


=  coskL-  TsinkL 

(All) 

The  complex  wave  propagation  constant  can  be  written 

k=  <pw2/E*)"’''2 

(A.  12) 

kL=uL(p/E'e'®)-'''2  =  “I-® 

(E'/p)'^^ 

(A.  13) 

(A.  14) 

=  wL  [  1  -  i  tan(-§)J 

<EVp)’''2sec(i)  ^ 

(A,  15) 

but. 

c=  the  phase  velocity  =  (EVp)’''^sec(-|) 

(A  16) 

C=  (a)L/c 

(A  17) 

Thus, 

kL=  411-  i  tan(-§))=  C  (1-iQ) 

2 

(A  18) 

where  O  =  tan(-) 
2 


(A  19) 


=  RkL 


7  = 


Mu' 


Mu' 


_  MkL 


AE*k  A(puVk2)k  Apt 

RCll-i?7l  where  R  is  the  mass  ratio 

=  M/pAL 


Plugging  in  7  and  k  into  (A.  1 1)  yields: 

coskL  -  TsinkL  =  cosC(l-iO)  -RC(I-IO)  sinC(l-IO) 
Expansions  for  the  complex  sine  and  cosine  are: 
cos(x-ly)  =  cosx  coshy  ♦  tsinx  sinhy 
sin(x-iy)  =  sinx  coshy  -  icosx  sinhy 


Using  the  above  identities  and  Eq.  (A.1 1) 
u(0)  Uo 

-  =  -  =  coskL  -  ?sinkL 

u(L)  u(L) 

=  cosCcoshCC7  ♦  isinCsinhCC?  -  RC(1-  iC7)  [  sinCcoshCO 
icosCsinhCO] 

This  can  be  expressed  in  terms  of  a  real  part, T^, and  an 
imaginary  part,  Tj ,  as: 

Tp  =  coshCO  (cosC  -  RCsinCl  ♦  RCO  cosC  sinhCf7 
and 
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(A.20) 

(A.21) 

(A.22) 

(A.23) 

(A.24) 


(A,25) 

(A.26) 


1’,  =  sinhCO  IsinC  ♦  RCcosCI  ♦  RCO  sinC  coshCO 


(A27) 
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